Strongly correlated dynamics in multichannel quantum RC circuits 
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We examine dissipation effects in a multichannel quantum RC circuit, comprising a cavity or 
single-electron box capacitively coupled to a gate and connected to a reservoir lead via several 
conducting channels. Depending on the engineering details of the quantum RC circuit, the num- 
ber of channels contributing to transport vary, as do the form of the interchannel couplings. For 
low- frequency AC transport, the charge-relaxation resistance (R q ) is a nontrivial function of the 
parameters of the system. However, in the vicinity of the charge degeneracy points and for weak 
tunneling, we find as a result of cross-mode mixing or channel asymmetry that R q becomes universal 
for a metallic cavity at low temperatures, and equals the unit of quantum resistance. To prove this 
universality we map the system to an effective one-channel Kondo model, and construct an analogy 
with the Coulomb gas. Next, we probe the opposite regime of near-perfect transmission using a 
bosonization approach. Focussing on the two-channel case, we study the effect of backscattering at 
the lead-dot interface, more specifically, the role of an asymmetry in the backscattering amplitudes, 
and make a connection with the weak tunneling regime near the charge degeneracy points. 

PACS numbers: 73.63.Kv,72.15.Qm,71.10.Ay 



I. INTRODUCTION 

The manipulation of mesoscopic systems to engineer 
quantum circuits has immense potential for future ap- 
plications. These unique systems exhibit a spectrum of 
novel phenomena which necessitates a better understand- 
ing of their dynamics. Technological advances have pro- 
vided the means to couple these systems to capacitive 
gates, thereby enabling detailed exploration of electronic 
transport at the nanoscale. 

In particular, the phenomenon of Coulomb blockade 
offers an excellent tool for the observation of interaction 
effects at the nanoscale and has become one of the 
cornerstones of modern condensed-matter physics. One 
of the simplest mesoscopic systems exhibiting Coulomb 
blockade is the single-electron box (quantum dot) and al- 
most two decades have elapsed since the first experimen- 
tal evidence of macroscopic charge quantization 4 . Quan- 
tum coherence and interaction effects drastically affect 
the properties of these systems. In fact, dissipation and 
resistance in single-electron boxes have recently sparked 
a growing attention in several parametric regimes both 
theoretically^ - — and experimentally^— . Specifically, 
the linear charge response to a gate voltage oscillation 
for single-electron devices has gained considerable atten- 
tion. 

In this paper, we investigate a multichannel version 
of the quantum Resistance-Capacitance (RC) circuit. A 
possible realization of this system involves a Coulomb- 
blockaded quantum dot coupled via a quantum point con- 
tact to a two-dimensional electron gas, which is in turn 
capacitively connected to a back gate (see Fig. []}. By 
tuning the opening of the quantum point contacts- 
through an auxiliary gate voltage, one can suitably 




FIG. 1: Schematic of a multichannel quantum RC circuit 
where M = 3 and M = 2, where M denotes the number 
of channels in the lead and Af denotes the number of chan- 
nels in the dot. All the possible lead-dot couplings are de- 
picted by black lines. The AC drive is capacitively coupled to 
the single-electron box. The geometric capacitance C g , is to 
be distinguished from the mesoscopic capacitance Co of the 
quantum RC circuit. 



control the number of conduction channels transmitted 
through the cavity. Here, we examine the low frequency 
behavior of the charge relaxation resistance in the cavity, 
in the limit of low temperatures such that quantum co- 
herence is preserved- For a small box, the role of coher- 
ence in charge quantization for a box consisting of non- 
interacting electrons has been studied 2 ^. In this work we 
investigate the role of Coulomb interactions in the charge 
quantization of the multichannel quantum RC circuit in 
the opposite regime of a large dot. For a single conduct- 
ing channel (mode), the relaxation of the charge on the 
dot when subject to an AC drive voltage has been studied 
by Biittiker et. al£-~—. The authors were the first to pre- 
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diet the universality of the charge relaxation resistance 
(R q ) in the context of a small (coherent) cavity, where 
the Coulomb interactions have been treated within the 
Hartree-Fock approximation. The quantum RC circuit 
has been realized experimentally in a two-dimensional 
electron gas, and the predicted charge relaxation resis- 
tance R q — h/(2e 2 ) has been confirme d 18 ' 19 ! 21 . 

Recently, the robustness of this value of quantized re- 
sistance in the presence of interactions in the cavity has 
been rigorously prove d 10 ' 11 . It has been shown that there 
is a mesoscopic crossover of the charge relaxation resis- 
tance from R q = h/(2e 2 ) to R q = h/e 2 as the size of 
the quantum dot is increased^ -. The value R q = h/e 2 
for a large (metallic) cavity in the vicinity of a charge 
degeneracy point has been obtained via an analogy with 
the Kondo mode l 29 ' 30 and is a consequence of the emer- 
gent Fermi liquid ground stat o 13 ' 31 " — . In this regime, an 
electron entering into the cavity is disentangled from an 
electron escaping the cavity. This unit of resistance can 
be viewed as two Sharvin-Imry contact resistances in se- 
ries. The charge relaxation resistance effectively captures 
interaction effects in a variety of exotic systems, such as 
fractional quantum Hall edge state o 10 ' 11 ' 34 , the Ander- 
son impurity mode l 11 ' 13 and topological insulator edge 
states 35 , all of which can be used to construct a quantum 
RC setup. The dynamical charge response in the case of 
a Majorana Coulomb box has also been analyze d 36 ' 37 . 

For a metallic cavity increasing the number of conduct- 
ing channels through the constriction causes the trans- 
port away from the charge degeneracy points to devi- 
ate from h/e 2 , and the value of R q is very sensitive to 
the engineering details of the system. The central mes- 
sage of this paper is that close to the charge degeneracy 
points however, when the tunneling couplings are weak, 
the charge relaxation resistance becomes universal and is 
in fact independent of the number of channels. To prove 
this result we project out states apart from those in the 
immediate vicinity of the charge-degeneracy point and re- 
formulate the problem in terms of a pseudospin^r— , and 
show via an analogy with the Coulomb ga a 43 ' 44 that the 
system can be mapped to the single-channel anisotropic 
Kondo model. 

The dynamics of the system in the opposite regime 
of perfect or near perfect transmission, i. e. strong tun- 
neling, can be studied using bosonization techniques. 
Adopting a generalization of well-known bosonization 
methodsi .^^ .^^ 7 . to Af channels, we show that R q 
is non-universal. For reflectionless channels with equal 
tunneling amplitudes R q — hj (A/"e 2 ). We focus on the 
special case Af = 2, and study the effect of backscat- 
tcring, in particular asymmetry in the reflection ampli- 
tudes, on the value of R q . For the Coulomb blockaded 
dot, the charge on the dot is pinned. We analyze the 
effect of charge fluctuations in the dot on the value of 
Rq for weak backscattering at the dot-lead interface. We 
show that a second order calculation in the backscatter- 
ing amplitudes, valid at high temperatures and frequen- 
cies, does not correct this value of R q . However, based 



on renormalization group arguments we conclude that 
the value h/e 2 of R q reemerges at low-energies, where 
the system flows to strong backscattering, in agreement 
with the weak tunneling analysis. More precisely, the 
system can be mapped onto an anisotropic two-channel 
Kondo model at the Emery-Kivelson line 4 ^,. An asym- 
metry in the backscattering amplitudes near the charge- 
degeneracy point causes the system to flow to a one- 
channel Kondo model at low-energies, since the channel 
with a stronger (bare) backscattering amplitude is even- 
tually perfectly reflected and pinched off 40 . 

In a different context, the variation of R q of a Coulomb 
box for a Landau-Zener sweep of the gate voltage has 
been studied by mapping the system to a dissipative par- 
ticle confined to a ring— li 5 -. For this out-of-equilibrium 
situation it has been shown that when a quantum of mag- 
netic flux is passed through the ring, the average of the 
relaxation resistance R q = h/e 2 . 

The remainder of the paper is organized as follows. 
In Sec. II, we present our results away from the charge 
degeneracy points via a perturbative expansion in the 
tunneling Hamiltonian following Ref. 10. In Sec. Ill, we 
discuss the underlying Kondo physics in the vicinity of a 
charge degeneracy point and formulate an analogy with 
the Coulomb ga a 43 ' 44 . In Sec. IV, we briefly address the 
situation at and close to perfect transmission. Appen- 
dices are devoted to technical details and mathematical 
derivations. 



II. WEAK TUNNELING ANALYSIS AWAY 
FROM THE CHARGE DEGENERACY POINTS 

The Hamiltonian of the system is given by 

H = ^ e k d\ a d ka +J2 £ p c \n c v& + E c(N - No) 2 

ka pj3 

+ w ( d i a c PP + c U dka ) ■ w 

kpafl 

Here, the subscripts a(= 1, . . . , Af) and /3(= 1, . . . ,A4) 
denote the channel index in the dot and lead respectively. 
At a general level, the channel (mode) index can also ac- 
count for spin degrees of freedom keeping in mind that 
tunneling between different spin projections is inadmis- 
sible. Hereafter, we consider the limit of large metallic 
cavity with a dense spectrum of energy levels. The sin- 
gle electron eigenfunctions of the two-dimensional system 
and the quantum dot are labeled by the index p and k 
respectively. The charging energy of the dot is expressed 
in terms of the charge operator 

Q = eN = e^dl a d ka (2) 

ka 

on the dot and Nq = C g V g /e is imposed by the gate volt- 
age (see Ref. llOi h We allow for a variety of inter-channel 
and intra-channel tunneling amplitudes (t a p), depending 
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on how the system is engineered. Note that in the follow- 
ing, we work at very low temperatures (T — > 0) in order 
to preserve the quantum coherence. 

In this Section, we derive the dependence of the relax- 
ation resistance R q on the number of channels (A4,JV) 
in the limit of weak tunneling amplitudes away from the 
charge degeneracy points of the dot, as given by Eq. ([lj. 
We assume that the charging energy E c is the most dom- 
inant energy scale, which implies the Coulomb blockade 
limits. The computation is a generalization of the scheme 
presented in Ref . [l(| • In this regime it is useful to group 
the terms in the Hamiltonian as follows, 

H = H' + H T , (3a) 

H = £ fc d{ a d ka +^2 £ p c lf}°pl3 + E c(N ~ ^o) 2 ' 

ka pa 

(3b) 

ht = w ( d i a c PP + c iv dka ) ■ ( 3c ) 

kpafi 

and treat the tunneling term Ht as a perturbation. 
The channels of the dot are denoted by a(= 1, . . . , Af), 
whereas for the lead it is given by j3(— 1, . . ., At). We 
study the AC response of the circuit using linear response 
theory. The gate voltage can be separated into its AC 
and DC components 



V g (t) = V g0 + V gl (t). 



(4) 



In the presence of a small time-dependent perturbation of 
the gate voltage, the charge on the dot Q — e(N) obeys, 



Q(w) = e 2 K(uj)V g (u), 



(5) 



where the retarded response function, following standard 
linear response theory 



K(t - t') = i0(t - t') ( [N(t),N(t') 



(6) 



measures the charge fluctuation induced by the AC com- 
ponent of the gate voltage. In the absence of electron tun- 
neling, the cavity charge in the ground state is (N) = N* 
and does not fluctuate, hence K = 0. We are at freedom 
to set N* = 0. Assuming weak tunneling, the charge 
fluctuations on the cavity are determined using pertur- 
bation theory in Ht- At T — the retarded and the 
time-ordered Green's function have a simple relation in 
frequency domain, so we instead compute 



K{t-t') = i(T t [N{t)N{t')] 



(7) 



We use the perturbative expansion in the interaction pic- 
ture 

K{t) = »<T t [JVj(f)JV>(0)]) 

_ . (0cs|Tt[JV7(t)#j(O)EMoo, -°*gs) 



{^Gs\T t [Ui{oo,-oo)]\(j) GS ) 



(8) 



Since the unperturbed Hamiltonian H' is not quadratic, 
Wick's theorem is not applicable. We therefore expand 
the evolution operator 

°° 1-00 pti 

£//(oo, — oo) = (— i) n I dt\ I dt2 ■ ■ ■ 

... r 1 dt n HT(h)H T (t 2 ) . . . H T (t n ), (9) 

J — OO 



in powers of H T (t) = e lH l H T (0)e 



-iH't 



In the ground 



state N\4>qs)=0. Thus, the zeroth and first order contri- 
bution to K (uS) vanish. The leading order contribution 
arises at second order. Using equations of motion the 
time-evolved electron annihilation operators can be ex- 
pressed as 



c P p{t) = e-^Cpp, (10) 
dka (t) = e-^+^^-^+^d^. (11) 

The charge relaxation resistance at low frequency is 
given by the expression: 



V g (oj) 



C (l + iojC R q ) + 0(oj 2 ), 



(12) 



the details of which are given in Appendix [X] The meso- 
scopic capacitance Co , which is distinct from the geomet- 
rical capacitance C g , in the weak-tunneling limit gives 



_ VQVlCg \ - 2 

C °~ 1/4- N§^ ta ^ 



(13) 



where i>q and vi represent the density of states in the 
lead and in the metallic cavity, respectively. Here we 
have retained terms to second-order in Ht- A fourth- 
order computation Ht allows us to extract the leading- 
order dissipative (purely imaginary) contribution to the 
function K(cj) in Eq. JS}. Thus, we obtain 



Ro 



2irh 

e 2 



(14) 



This illustrates that in general, away from the charge 
degeneracy points, i.e., for Nq ^ 1/2, the charge relax- 
ation resistance R q is not universal as it depends on the 
engineering details of the system. For the case of diago- 
nal couplings tap = t a 5 a fj (J\f = M), R q reduces to the 
non-universal result 



Ra 



h 
72 



V t 4 



(15) 



It should be noted that in the case of (almost) isotropic 
diagonal couplings, i.e., t a = t, the quantum RC cir- 
cuit becomes dissipationless in the limit Af — > oo as 
R q = h/(e 2 AT). This value demonstrates the violation 
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of Kirchhoff's law, which in contrast predicts that R q 
would be inversely proportional to the sum of the trans- 
mission probabilitie s 18 ' 19 ; and is a direct consequence of 
the quantum (phase) coherence. 

We remark that for the special case when the trans- 
mission prabability of a single channel dominates ,i.e. 
ti 3> {t2, ■ ■ ■ , tj\f} one recovers the value R q = h/e 2 . 



III. KONDO MODEL AND COULOMB GAS 
NEAR THE CHARGE DEGENERACY POINTS 

The weak-tunneling analysis of Section II leads to the 
surprising conclusion that for the case of large inter- 
channel mixing R q = h/e 2 , i.e. a reemergence of the 
unit of resistance. This follows directly from Eq. (fT"4")) 
if we take the amplitudes of the tunneling matrix to be 
approximately of equal strength, i.e. t a p f» t. It is possi- 
ble to reformulate the quasiparticles of the dot and lead 
(individually) in terms of the totally symmetric combi- 
nation of the modes and additionally construct TV — 1 
and AA — 1 mutually orthonormal modes respectively. In 
this new formulation t a p = t is the condition for perfect 
reflection of the additional TV— 1 and AA — 1 modes respec- 
tively, and effectively only a single channel (the totally 
symmetric mode) is transmitted. The value R q = h/e 2 
is a consequence of this emergent one-channel quantum 
RC circuit. 

Due to the strong Coulomb blockade in the vicinity of 
a charge degeneracy point TVo = n + 1/2 (n 6 Z), we 
can project out all charge states other than those corre- 
sponding to n and n + 1, which then mimic a spin- 1/2 
particle 3 ^. We show that this gives rise to an emergent 
one-channel Kondo model and also arrive at this result 
by mapping the system to the Coulomb ga n 43 ' 44 . 



A. Mapping on a Kondo model 

For the Coulomb blockaded dot near the charge degen- 
eracy points No = n + 1/2, the effective Hamiltonian of 
the system can be written as^ 8 . 

H = ^ e kd\ a d ka + ^ £ p c I(s c pI3 ~ h z S z 

ka pj3 

+ W ( d la c Pf> S ~ + cl? dkaS+ ) ■ ( 16 ) 

kpafi 

The effective magnetic field h z = e [Nq — (n + 1/2)] /C g 
denotes the deviation from the charge degeneracy 
points^ 8 - and is assumed to be small compared to E c = 
e 2 /2C g . Let us consider purely diagonal couplings, i.e., 
for t a p = t a 6 a /3 (Af = AA.). It is straightforward to 
see from Eq. (|16p that one recovers the anisotropic Af- 
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unrestricted 




4-channel mixing 




3-channel mixing 




— 2-channel mixing 




no channel mixing 





2 4 6 8 10 12 

In (E c /T) 

FIG. 2: Renormalization of the DC conductance for a multi- 
channel quantum RC circuit (M = Af = 5) with decreasing 
temperature, at the charge-degeneracy point. We consider 
randomly generated couplings (unrestricted) and eliminate 
progressively the interchannel-couplings, starting with the 
farthest channels(4-channel mixing) till all the inter-channel 
couplings are removed(no channel mixing). The random cou- 
plings are normalized according to the degree of channel mix- 
ing, such that g(0) = 0.5. The dashed curve corresponds to 
the unitary limit of the Kondo model. 

channel Kondo model 

H = ^ ^kdkadka + £ P C p/3 C P0 ~ hzSz 
ka pP 

+ ^2 ta ( d ka C pa S ~ + 4>a dkaS+ ^) ' ( 17 ) 
kpot 

Introducing a pseudospin index (er = J [, \.) which denotes 
the position of a particle in the system, reservoir lead (f) 
versus cavity (J,). Furthermore, defining dka — o-k^a and 
Cpa = o-kia we can rewrite the Hamiltonian of the system 
as 

-Hdiag = ^2 / Ek IJ a\ iya akaa ~ h z S z 

kaa 

+ Y, ta ( a ltc* a kl<xS~ + a{i a a k^S + ^j , (18) 

kpa 

which is precisely a TV-channel Kondo model with J± a — 
t a and J|| a = O 4 ^— . Here, we assume that — tk — 
Eki- A small deviation from the resonance condition 
plays the role of a magnetic field, leading to a Zeeman 
term in the Kondo model. 

This equivalence with the TV-channel Kondo model 
leads to the value R q = h/(Afe 2 ) and corroborates the 
prediction of Eq. (fT5)) in the weak-tunneling limit away 
from the charge-degeneracy points. The two-channel sit- 
uation has been studied in detail by two of us 3 ^. 

As emphasized in Ref. Hl|, the general situation in 
the single-electron box is fairly intricate. In particular, 
cross-mode (channel) tunneling is a relevant perturba- 
tion which leads the system ultimately to a one-channel 
Kondo fixed point in the case of spin-polarized electrons. 
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This follows from the low-temperature effective Hamilto- 
nian given by Eq. (| , which predicts the emergence of a 
unique effective tunneling mode in the lead and the elec- 
tron box. The other modes are perfectly back-scattered 
when the inter and intra-channel couplings are identical. 
Deviations from this condition causes partial transmis- 
sion of these other M. — 1 and N— 1 channels respectively, 
and leads to corrections of R q . However, renormalization 
group analyses show that these perturbations are irrele- 
vant, and for low energies the value of R q is stable^. The 
energy scale A* which determines the crossover from Af- 
channel to 1-channel regime, depends on the precise form 
of the tunneling amplitudes. When the off-diagonal ele- 
ments are small and/or sparse, the onset of the effective 
one-channel behavior is further delayed, i.e., one has to 
go to lower temperatures to measure R q — h/e 2 . 

In Fig. [5] we show the variation of the DC conductivity 
of the multichannel quantum RC circuit, as a function of 
temperature. The conductivity G(T) = 47r /t e g, where 
g (hi(E c /T)) is a dimensionless parameter, has been nu- 
merically obtained using the scaling equations in Ref. 
l4ll . We fix the number of channels Af = Ad = 5 in 
Fig. [21 and consider a randomly generated set of tunnel 
couplings (magenta curve). Using this matrix of tun- 
neling amplitudes, we then restrict ourselves to only the 
diagonal values (blue curve), setting all the off-diagonal 
entries to zero. The tunneling amplitudes are globally 
scaled such that g(0) = 0.5, corresponding to T = E c . 
Following this procedure, we then allow off-diagonal di- 
agonal elements representing mixing of adjacent chan- 
nels (green curve). We then include the elements cor- 
responding to mixing of adjacent and next-to-adjacent 
channels (red curve), and ultimately additionally include 
the terms represented cross-mode tunneling of next-to- 
next-to-adjacent channels (gray curve). From Fig. [21 it 
can be seen that the system behaves identically to an ef- 
fective one-channel quantum RC circuit, at low tempera- 
tures. Asymmetry in tunneling amplitudes and/or inter- 
channel mixing are the factors responsible for this emer- 
gent one-channel behavior. As apparent from the curves, 
introducing greater inter-mode couplings, increases the 
energy scale A* (given by the minima of the curves) at 
which the system crosses over to the one-channel regime. 
The (black) dashed lines correspond to the unitary limit 
of the 1-channel Kondo model. The effective Kondo tem- 
perature Tj^ 5 - , is the point at which the curve converges 
to the unitary limit. For purely diagonal and equal cou- 
plings, the system is equivalent to the A/"-channel Kondo 
result. The crossover energy scale A* as a function of 
the number of channels, has been studied 4 ^, and shown 
to be suppressed as the number of channels is increased. 

The one-channel Kondo fixed point is stable whereas 
the A/"-channel is unstable, and a deviation from the con- 
dition t a p = tS a ff causes the system to flow away from 
the A/"-channel fixed point. This implies that the quan- 
tized resistance unit R q = h/e 2 may be observable in 
a multichannel quantum RC circuit, if it is possible to 
access the energy regime below A* experimentally. 



Deviations from the charge-degeneracy point, i.e. an 
effective magnetic field in the pseudospin language, may 
inhibit the onset of the effective one-channel behavior. It 
is a relevant perturbation and the renormalization flow 
needs to be stopped when the cutoff (temperature) equals 
the effective magnetic field, a condition determined self- 
consistently 41 . If this happens before the crossover to the 
one-channel regime, then the universality of R q is killed. 
The variation of R q for the case of an infinitesimal field, 
such that the energy scale lies between Tk and A*, has 
been studied perturbatively and the one-channel result 
has been confirmed 56 . However, the value of R q precisely 
at the step, i.e. below Tk, remains an open question. 



B. Connection with the Anderson- Yuval expansion 

In the previous Subsection we showed that below the 
energy scale A*, the multichannel quantum RC circuit 
near a charge-degeneracy point is equivalent to a one- 
channel Kondo model. 

Here, we provide a transparent way to better under- 
stand the connection of the quantum RC circuit with the 
one-channel Kondo model, based on an analogy with the 
Coulomb ga a 43 ' 44 . We assume that a strong Coulomb 
blockade exists in the quantum dot, and for the sake 
of simplicity restrict ourselves to the charge-degeneracy 
point. Deviations from this point, assuming that the Zee- 
man energy is negligible compared to E c , can be readily 
incorporated. However, this makes the mathematical de- 
tails significantly more complex without being more il- 
luminating, and will hence be ignored for the sake of 
clarity. Using the Hamiltonian in Eq. ([T]) we can write 
the partition function as the imaginary-time functional 
integral 



[Dd][Vc]e 



(19) 



where the action 



S 



dr 



S^dkc 

ka 



[d T + e k }d ka +E c (N-N ) 2 + 



^ C p p[d T + S p ]c p p + t ^2 (cppdka + d"kaCp/3 



(20) 

p/3 kpa/3 

In the above action, we have assumed that we are below 
the energy scale A*. In this regime, it is possible to select 
an appropriate unitary transformation of the degrees of 
freedom of the system {c ka , dpp} — > {c ka , d p ^}, such that 
the channel- mixing is maximal (t a p = t) in terms of these 
effective degrees of freedom. 

Introducing the Hubbard-Stratonovich variable V(r) 



via the identity J[DV] exp - J f 



, ^dr 



1, the quar- 



tic term in the d kc 
transformation 



operators can be absorbed by the 



V(t)^V(t) + 2iE c (N-N ), 



(21) 
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such that 



[Vd][Vc\[VV]e- 



(22) 



Here, the transformed action 



S 



(It 



d ka [d T + e k + iV]d ka + — - iiV T 

ka C 



+ X! CpPftr + £ p\C'PP + * X! (ppadk(3 + h.C 
p/3 kpa/3 



(23) 



The variable V(r) can be expanded in terms of its Fourier 
components V(r) = ^ ^K n e"'' mT , where 2/ m = 22p de- 
notes bosonic Matsubara frequencies. The zero mode 
Vo can be pinned to the values 2nn, where n G Z, 
and its fluctuations neglected^,. We define the field 
<f>n{T) = fn V n (T')dT' such that it satisfies the bound- 
ary condition </>„(r + /3) = 4> n ( T ) + 27m. Transforming, 



dk a (r) -t d fcQ (r)exp [-i(/)„(r)] 



(24) 



we decouple iV(r) and V(r), while preserving the an- 
tiperiodicity dfc Q (r + /?) = — dfc a (r). Thus, we obtain 



oo . 

Z =H \ [Vd)[Vc][V<t> n ]e- s , 



(25) 



where the transformed action is given by 



S= dr 
Jo 



^2 dkaldr + £k]dk a + —{dr^nf - iN d r 

ka ° 

+ X] ^ + ^ fc 
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(26) 



Integrating the topological term Nod T (f> over (imaginary) 
time we obtain 



[I'd] [2?c] [£></>„] exp[— S] . (27) 



We separate the terms in the action as S = Sq 1 +Sq +Si nt , 
where we define 



Sf= I dT 



(28) 



the action for noninteracting electrons in the lead and 
dot, the kinetic energy of the field <f> 



1 

S ° = I dT AE r V dr J ' 



(29) 



and finally the interaction term 



dr t V ( 



^ K Cpffd ka e-^ +h.c.) . (30) 

/cpa/3 



Expanding the partition function in powers of Si n t, and 
subsequently transforming to the charge representation 
we integrate out the fermionic degrees of freedom. The 
existence of a strong Coulomb blockade together with 
the fact that we confine ourselves to the vicinity of 
the charge-degeneracy point, allows us to project out 
the other charge states which are energetically distant. 
This leads to the partition function of the well-known 
Coulomb gas^ii 



Z = J2 



wit 2 ATM) 



1=0 





r 1 -^ d T2 


j 




Jo T c 


Jo 



dT, 



21 



■ exp 



2 log 



(31) 



i<j 



Details of this procedure are included in Appendix |B| 
Here, the ultraviolet cutoff l/r c is formally identified to 
be A*. 

To make a direct connection with the Kondo mapping 
in the previous Subsection, we could have alternatively 
rewritten Eq. (|20j) in terms of the totally symmetric mode 
and Af — 1 and M. — 1 modes in the dot and lead respec- 
tively. The charging energy term is given exclusively in 
terms of the totally symmetric mode, and the Coulomb 
repulsion term is what enforces the projection to the 
two charge states near the charge-degeneracy point. The 



Coulomb gas expansion follows directly from this fact. 
This Coulomb gas representation allows us to prove that 
at energy scales smaller than E c , the effective model for a 
quantum RC circuit with maximal inter-channel mixing, 
is the one-channel Kondo model. 

In the complete absence of channel mixing, i.e. diag- 
onal tunneling matrix elements, rewriting the partition 
function as a cumulant expansion in Ht and retaining 
terms to second-order yields the well-known effective ac- 
tion for the single-electron box in terms of a dissipative 
particle on a ring 58-64 . The key ingredient in this map- 
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ping is the fact that higher orders in the cumulant ex- 
pansion are suppressed by factors of at least 0(1/ AT). 

For large- Af this furnishes a natural hierarchy of rele- 
vant terms, and in the Af oo regime all other terms 
become completely irrelevant. This is in contrast to the 
case where channel mixing or anisotropy in the tunneling 
amplitudes is present, where no such ordering is possible. 
Defining k — (yQV\t 2 AfAA), all terms in the cumulant 
expansion in k are of 0(1) in Af and A4, once k is held 
fixed. 



The kinetic term takes the form 40,46 : 



v F 



#k = ^ / dxJ2 [{d^i{x)f + {d x 6i{x)f]. (32) 

Zn J -°° 1=1,2 

Here, vp is the Fermi velocity which is obtained by lin- 
earizing the energy spectrum around the Fermi points. 
The ultraviolet cutoff 1/et denotes the extent around kp 
upto which the spectrum can be linearized. We note that 
the bosonic field operators (x) and 9i (x) obey the usual 
commutation relations: 



-sgn(x - y)Si 



(33) 



IV. CLOSE TO PERFECT TRANSMISSION 

We now address the regime opposite to that in Section 
III, i.e. a Coulomb-blockaded dot near perfect transmis- 
sion. In Sub. IIV Al we study the perfect transmission 
limit, where the system can be exactly solved. We intro- 
duce the notation and find the general expression for R q 
and Co for Af channels. We then focus (Sub. IIV Bp on the 
case of 2-channels, where we use an intuitive description 
of the system in terms of spin and charge degrees of free- 
dom, and introduce weak-backscatterin g at the dot-lead 
interface. We extend the results of Ref. |36|, which treats 
the case of perfectly symmetric channels, to include an 
asymmetry in the backscattering amplitudes. Such an 
anisotropy in the backscattering amplitudes for a single 
(orbital) channel of spinful fermions can be introduced 
by applying an in-plane magnetic field^. At high en- 
ergies we treat perturbatively renormalization effects of 
the backscattering amplitudes and show that the value of 
R q remains unaltered. Lowering the energy cutoff leads 
to non-universal corrections and we conclude that in the 
low-energy limit one of the channels is pinched off and the 
system flows to a single channel Kondo model 40 . Here, 
Rq equals the unit of quantum resistance, and near the 
charge-degeneracy point the physics is identical to the 
one-channel Kondo model obtained in the weak-tunneling 
regime discussed in Section HT1 



A. Perfectly transmitting channels 

It is convenient to describe the system using a 
bosonization approach following Ref. I4al46l . The coor- 
dinates of this one-dimensional system are chosen such 
that the lead occupies the semi-infinite line x £ (— oo, 0). 
The dot extends from x — to x = L. We consider the 
situation when the level spacing of the dot A = (ttvf)/L 
goes to zero, which implies that L — > oo. In this limit 
electrons entering and leaving the dot are uncorrelated. 
For a concrete analysis, we focus on the case with two 
conducting channels. We bosonize the total Hamiltonian 
H = H k + H c + H hs m (for an excellent survey on the 
subject and details of the notation used, see Ref. l39h ■ 
Here, we assume that the channels are perfectly tansmit- 
ting. The role of backscattering will be investigated in 
the next subsection. 



In particular, the total electron charge on the cavity be- 
comes Q = (l/7r)(</>i(0) + 02(0)) such that the charging 




&(Q)-7riVb 



(34) 



and No, which is related to the gate voltage, has been 
introduced earlier in Sec. II. For electrons in the cavity, 
we choose </>i(oo) = which fixes the charge on the dot to 
be zero when Nq = 0. Here, the two conducting modes 
can refer to the two spin polarizations of the electron or 
two electron channels. In either case it is then natural 
to introduce the charge and (pseudo-)spin modes cf> C)S = 

(M0)±M0))/V2.^ 

At perfect transmission, we can omit the spin part of 
the Hamiltonian and integrate out the charge fields ex- 
cept the one at i = 0, which results in the following 
action: 



s c = ^U c (o,w„) 



^ ) • 05) 



Here, we have defined the shifted charge field 4> c (0, w„) = 
</> c (0, uj n ) — ttNq, in order to make the action independent 
of N . From the Green's function of the charge mode 
4> c (fi, u n ), one can easily evaluate the function K(t) in 
Eq. ([7]) and obtain 



( 1 



C' 



1 _ IUJTX 

2E C 



(36) 



Above, we have chosen units in which H = 1. Comparing 
the low-frequency expansion of Eq. ([31]]) with Eq. (TT2")) 
this results in: 

h 



Rq ~ 2e 2 ' 



(37) 



At a general level, for a large cavity and a number Af of 
perfectly conducting channels one finds R q = h/(Afe 2 ). 
In Section III we showed that the presence of channel 
mixing and/or asymmetry in the tunneling amplitudes 
cause Rq to be universal and equals the unit of resistance 
h/e 2 . In the remainder of this section we study the effect 
of backscattering on R q ; in particular, we investigate the 
role of asymmetry of the backscattering amplitudes on 
the relaxation resistance. 



s 



B. Weak backscattering and 2-channel Kondo 
physics 

A deviation from the perfect transmission condition is 
captured by the back-scattering Hamiltonian 40 



H bs = — V | A, | cos (20,(0)). 



(38) 



In terms of the charge and spin fields this can be recast 
as 

H hs = —y"\\ l \cos(V2{(f> c + Vl <f> s )), (39) 



where we have introduced r\\ — 1 and r\2 = — 1. Defining 
A = | Ai| + j A2 1 and 5\ = | A2 1 — |Ai| we can rewrite Eq. 

as 



A cos fV2^ c — nNnj cos ^\/2</> s ^ 
+ S x sin U/2<j> c - nN^j sin \y/2<j>^ . (40) 



In the Coulomb blockade regime, the shifted field <\> c is 
pinned to to zero. To be mathematically rigorous we 
adopt a coherent state functional integral approach in 
the ensuing treatment. The reader is to assume that the 
various transformations are implemented in the action 
and not at the level of operators. 

' A cos(7rJV ) \ ' 
/2-na J 



We define the energy scales F_| 



VF \ 



and T_ 



VF V 



<5 A sin(nN ) 



A physical interpretation of 

these energies in terms of lifetimes of Majorana fermions 
is discussed later in this subsection. To study the low- 
energy physics, following Ref. [3II we introduce an inter- 
mediate energy scale A, such that max{r+, L} « A < 
E c . We then integrate out the high-energy contribution 
A < to < E c to the charge mode, such that the resultant 
action still effectively describe the interaction between 
charge and spin modes at the energy scales of T± . This 

amounts to replacing the factor cos ( \/2(f> c — ttNq ) with 



cos (V24> l c - ir N Q 



(41) 



where <jr c contains the low-energy fluctuations of the 
charge field. In this regime, we can expand the backscat- 
tering terms in powers of (fr c , and approximate 



cos 



W2$ c - ttNqJ « cos(ttAo) + V2sia(nN a )$ l c 

- cos(7rA )(^) 2 := cos(ttAo) + V2mC + (42) 
in (V2$ l c - ttA ) w - sin(7rA ) + V2 cos(ir N )<f>' c 
+ sin(7rA )(^) 2 := -sin(7rA ) + V2^C_. 



We define the dimensionless backscattering amplitudes 
n 2 via the relation 



A1.2 = ri 



2avpEc^ 



(44) 



:( T )0c(r)) £ >A 

vp/ja, 7 being 



In other words we set 
(1/2) ln(uj£)7r/2E c ), where ujjj 
the Euler-Mascheroni constant. 

In the ensuing analysis we implicitly assume (j> c to be 
small. Since the charging energy E c is assumed to be 
much larger than the backscattering amplitudes, the dy- 
namics of the spin field are completely captured by the 
pinned charge field <jr c . Next, we compute corrections to 
R q due to charge fluctuations and study the role of an 
asymmetry in the backscattering amplitides. 

The spin part of the Hamiltonian can be treated ex- 
actly by refermionizing the field 4> s using the procedure 
outlined in Ref. H(J The kinetic part of the Hamilto- 
nian can be written as the sum of its charge and spin 
parts, i.e. i/k = + H^. Rewriting H£ as a functional 
integral and integrating out the charge field away from 
i = 0we obtain the action in Eq. (|35|) . In terms of the 
new fermionic fields %p s [x) and the "slave" fermion d, the 
kinetic part of the spin Hamiltonian becomes 



iv F / dxtpl(x)d x 'ilj s (x), 



(45) 



and the backscattering Hamiltonian is given by 



(A+ + A(7 + )(v4(0) + ^(0))(d-d+) 

f (x_ + s x 6-) (v4(o) - ^,(0)) (d + S) 



(46) 



reminiscent of the two-channel Kondo mode l 48 ' 67 in the 
presence of a channel asymmetr y 40 ' 68 . The introduction 
of two Majorana fermions follows from Ref. 0, and Ref. 
68 uses this representation to make an analogy with the 2- 
channel Kondo model in the presence of a channel asym- 
metry (at the Emery-Kivelson line). This Majorana model 
also finds applications in the context of dissipative meso- 
scopic structures^. 

Here, we have defined X + = \ cos^ N$) / and 
A_ = —8\ sin(7rAo)/v / 27ra, to isolate the terms which in- 
volve the shifted charge field <jr c . These terms capture the 
dependence of the spin propagator on the pinned charged 
and is included in the unperturbed Hamiltonian. Next, 
the action of the remainder of H bs on the charge response 
function is treated perturbatively. For this purpose it is 
convenient to recast this Hamiltonian in terms of Majo- 
rana fermions on the lines of Ref. H3- Details of this 
computation are provided in Appendix [Cl 

We define the functions k + (No) := cos(ttNq) and 
K-(Nq) := sin(7rAo). The energy scale of the Fermi ve- 
locity, which is lost in the bosonization process where 

an infinite bandwidth is assumed, is reintroduced by fix- 
(43) 

ing the cutoff 45 . The widths T + and T_ of the impurity 
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2X 2 /v F , 



determine the two 



Majoranas, given by — ^i^/ 
intrinsic energy scales in the problem, corresponding to 
the total backscattering strenth and the channel asymme- 
try respectively; they can be explicitly obtained through 
the Renormalization Group (RG) approach^. We define 
the infrared cutoff V — max{r+, It corresponds to 

the energy scale which is inversely proportional to the 
lifetime of the (shorter lived) Majorana fermion. As pre- 
viously stated, the system can be mapped onto the 2- 



channcl Kondo model at the Emery-Kivelson linei^ and 
in terms of this they represent the Kondo coupling con- 
stants J x and J y respectively. 

The expression for the charge reponse function to lead- 
ing order in the backscattering can be written as 

K(u) = K q {uj) + Ki{uj) + K 2 (u) + 0(A 4 ), (47) 

where 



J 



K 2 (u) 



2nE r ( 

1 l 2E c 



J2 K d N o) 2 (ri + (r 2 f 



C=± 



log I — 



2wE r ( _ ^ 1 

1 *2E C 



5>_ c (7Vo) 2 (r 1+ Cr 2 ) 



1 + ,- log( 



r 



(48) 



When the backscattering amplitudes are symmetric the 
contributions corresponding to £ = — disa ppe ar and this 
case has been looked at in detail in Ref. 36. The level 
asymmetry, which is captured by the energy scale T_ 
plays a crucial role in the low energy behavior of the 
system which we discuss below. 
We now compare 



K{u)=i — 



dte 



" ({k(r)M0)}) (49) 



with the expansion K(uS) — Co(l + iujCoR q ) + (cj 2 ) . To 
proceed systematically we formally expand the functions 
Kq(uj), K\(uj) and K 2 (uj) to linear order in the frequen- 
cies 



Kq{ui) = A + iujB 



iujB[ 



(50) 



C=± 

From this, we find the mesoscopic capacitance 
C Q = A Q + ]T r 2 [A{ + Ai 



R a 



(A> + £ 



r 2 



A[+Ai 



(51) 



where we have retained terms to leading order in the 
backscattering amplitudes. Here, we have defined the 
parameters r± = r\ ± r-i. 



One then finds to 0(r 2 ) 



+ 



2ttE, 



(r 2 cos(2ttN ) log 



R « = 2^- 



(52) 



In the vicinity of No = 1/2 the logarithm term dom- 
inates and in particular setting r\ — r 2 one recov- 
ers the correction to the mesoscopic capacitance 6C = 
2-f/{-KE c )r 2 cos(27r7V ) log (l/(r 2 cos 2 (7rJV ))) obtained in 
Ref. |38|. Specifically, the logarithmic contribution of 
the capacitance Co agrees with Ref. |4CJ in the presence 
of a channel asymmetry. Note, R q is not affected by 
backscattering and remains h/2e 2 to the second order 
in the backscattering amplitudes. Non-universal correc- 
tions to R q nevertheless appear at fourth or higher order 
in backscattering amplitudes. 

At low frequencies and temperatures the backscatter- 
ing amplitudes are substantially renormalized, thus en- 
tering the regime of strong backscattering, and this per- 
turbative scheme to compute R q is inadequate. For low 
energies, the channel asymmetry for energy scales be- 
low r_ grows till r max = max{ri,r2} — > 1, indicating 
that one of the channels is pinched off and there is only 
one channel effectively contributing to transport^. In 
this regime, near a charge degeneracy point, the system 
is described by a 1-channel Kondo model and the value 
h/e 2 of R q should reemerge in agreement with the weak- 
tunneling limit. 
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V. CONCLUSION 

To summarize, we have investigated dissipation effects 
in a strongly correlated multichannel quantum RC circuit 
through the charge dynamics at low frequency and the 
concept of charge relaxation resistance. We have corrob- 
orated the violation of the Kirchhoff 's law as a result of 
quantum coherence effects. Our main message is that in 
the vicinity of a charge degeneracy point and in the case 
of a large cavity, the system can flow to a one-channel 
Kondo fixed point at low energy as a result of channel 
mixing or channel asymmetry, leading to a resistance unit 
h/e 2 close to absolute zero reminiscent of the one-channel 
RC circuit 1 ^. 

It is relevant to note that such a resistance unit has 
also appeared in a distinct purely non-equilibrium con- 
text of Landau-Zener sweep on the gate voltag o 14 ' 15 . In 
the latter situation, the effect of channel mixing or chan- 
nel asymmetry was not considered which suggests a dif- 
ferent origin for the emergence of the resistance quantum. 

It is also important to underline that for the multi- 
channel quantum RC circuit, prominent deviations from 
the resistance quantum are expected either in regions 
where the weak backscattering or weak tunneling limit 
describing the coupling between lead and cavity is ap- 
plicable. We believe that our results shed new light on 
the understanding of the dynamics of quantum RC cir- 
cuits in close connection with multichannel quantum im- 
purity model a 49 ! 70 , and could be in principle detected 
with current technolog y 2 ^ 30 , especially close to perfect 
transmission^!. 

We thank B. Horovitz and P. Le Doussal for fruitful 
discussions. This work is supported by Department of 
Energy under the grant DE-FG02-08ER46541 (P.D and 
K.L.H). 



Appendix A: Perturbation Theory in weak 
tunneling regime 

In this Appendix, we present the derivation related to 
the main Eq. (fT4")) in Sec. II. We start from the function 
Kit) defined in Eq. ([7]) and first we focus on the second 
order contribution in the tunneling amplitude, which will 
generate the mesoscopic capacitance Ccr&. We generalize 
the computation of Ref. 10 which has been performed in 
the context of the one-channel quantum RC circuit. 



1. 2nd order 

Assuming t > the response function can be simplified 
to give 

oo rO 



K (2 \t) =-i dh / dt 



x (<t>Gs\HT(h)N(t)N(0)H T {t 2 )\<j> GS ). (Al) 



The contributions can be conveniently indexed by the 
path adopted by the dot electrons. There are two possible 
paths which give non-zero contributions. — > 1 — > 
represents a lead electron hopping onto the dot and back, 
while — > — 1 — > represents a dot electron hopping onto 
the lead and back. 

For t > the contribution of the — > 1 — > path is 
evaluated to be 



dh 



kpa/3 



x 9(-e p )8(e k ). 
Similarly, the — >• — 1 — >• path contributes 



(A2) 



Kpap 



dti 



fa e -i[e p -e fc +B c (X+2JV )](ii-t2) 



x 9(e p )9(-e k ). 



(A3) 



Collecting these contributions together, we obtain for 
(OO) 



kpcxfi 



6(e k )0(-e p ) iAl (k,p)t 

e(-€ k )o(s p ) iA_i(fc, P )f 



(A4) 

Here, we have defined 

A ±1 (k,p) = E c (lT2N )±(e k -e p ). (A5) 



The t < contribution can be deduced from Eq. (| A4I) by 
the replacement t — > —t. Fourier transforming we obtain 

/•OO 

K<- 2 \uj) = / dte iut e-° + ^K^(t) 



t 2 fl(-g p )fl(efc) 



1 



A a (k,p) - lo + ir] 



(A6) 



Assuming a metallic dot, we can let J^ fc — > vi de k 

(for the lead we have ^2 p — > vq de p ) and the meso- 
scopic capacitance is determined to be 



Co = e 2 ^)(0) = T ^%5:^. 



1/4 -TV 2 



(A7) 



a/3 



This corresponds to Eq. (|T5|) in the main text. This 
shows that in the weak-tunneling regime, the meso- 
scopic capacitance of a single-electron box differs very 
strongly from the geometrical capacitance as a result of 
the Coulomb blockade phenomenon. 
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2. 4th order 

The various contributions which arise at 4th order are 
classified by indexing the dynamics of the charge occu- 
pation on the dot and the number of virtual particle-hole 
excitations. Only the paths 0— >1— >1— > 0, 
O^l^O^ -1 -> 0, -)• -1 ^ 1 ^ 
and 0^— 1^0— >— 1— >0 contribute at low en- 
ergies (i.e. linear in uj). Furthermore, based on phase 
space arguments, only paths within this subset that 
have particle-hole virtual excitations contribute. We ex- 
plicitly compute the relevant contribution by the path 
->■ 1 ->• ->• 1 ->■ 0. Its contribution to for t > 
follows from the expression 



ft rt 2 pO 

dh I dt 2 / dt 3 / dt4{^as\HT(ti)N(t)H T (t2) 



The contribution of this term to the response function 
K(w) to linear order in ui is evaluated to be 

<* 2 & 2 

9(e k ) 



[E C {1 - 2N ) + e k } 2 



+ 0(uj 2 ). (A13) 



Note that the imaginary contribution for uj < arises 
from the t < expression. Collecting the contribution 
from the paths O^l^O^ -1 -> 0, 0^-1^0^ 
1 — > and ^ — 1 — > — > — 1 — > in a similar fashion, 
and noting that for electron-hole excitations in the lead 
given by the grouping in Eq. (|A10|) we obtain an identical 
contribution 



t JO JO J-oo 

x H T (t 3 )N(0)H T (U)\<j> GS ) 



(A8) 



Isolating the product of operators contributing to this 
path we obtain 

poo pt pt 2 pO 

^ taxh^a^i^az^ot^i \ dt\ / dt 2 / dt 3 / dt 4 

„,„,«, Jt JO JO J-oo 



^smK^(uj) = 2ttUqV^uj ^ t ai p 1 t ct2 p 1 t a2 p 2 t ail 3 2 

a 2 3 2 

1 1 



E C (1-2N ) E C (1 + 2N ) 



0(u A ). 



(A14) 



...yS 4 



x e 



-i(Ai(fei,pi)ti-Ai(fc2,P2)<2+Ai(fe3,p3)t3-Ai(fc4,p4)t 4 ) 



X (^Gsl4 lft 4 lQl 4 2 a 2 c P2/32Cp3/33 dfc 3«34 4a4 c P4/34l? !, Gs)- 

(A9) 

The pairings corresponding to single particle-hole excita- 
tions in the intermediate state are given by 



In the last step we assume a metallic dot and let — ¥ 
Vl f X oo dek ( for tlie lead we nave S P v o J^dep). 
Comparing Eq. (|A14[) with the general relation — 

C (l + iuiC Q R q ) + 0(ui 2 ) and using Eq. ([AT|) we identify 
the AC charge relaxation resistance to be 



Ra 



2irh 

e "2 



{4>Gs\cl lPl d kiai d\ 2a2 c P2 p 2 c\ 3 03 d k3 Q3 4 4Q4 c P4/34 \4>gs) 



Pi tct 1 /3 1 ta 2 fl 1 ta 2 f3 2 taip 2 

<*202 



(A15) 



(A10) 



and 



{4>GS \Cp l/3l dk iai d\ 2Ct2 Cp 2 p 2 Cp 3/ 3 3 dk 3 a 3 4 4 a 4 C PiPi \4>Gs) 



(All) 



This corresponds to the formula (|T4l in the main text. 
Note that here we have explicitly restored the Planck 
constant h which was set to one in the derivation. 



Appendix B: Coulomb gas formulation 



The pairing (jAlOp corresponds to an electron-hole exci- 
tation in the reservoir, whereas (|A1 1|) corresponds to an 
electron-hole excitation in the dot. Integrating over the 
time variables and extracting the imaginary contribution 
of the above expression leads to a contribution linear in 
u> (hence gapless) 

tcnpJ>a.2P\ta.2pJ>a l p 2 Q (e fel ) 6 (e k2 ) 9 (— e pi ) 9 (e P2 ) 

kipik 2 p 2 
oi/3lQ 2 /32 

pi[Ai(fci,p 2 )-Ai(fci,pi)]t 

x (A12) 

Ai(*i,pi)Ai(fc2,pi)Ai(fca,j>a)Ai(fci,j)2) 



Here, we show that for energy scales below A*, the 
partition function close to the charge degeneracy points 
can be rewritten exactly as the partition function of the 
one-channel Kondo model, resorting to the Coulomb gas 
analog y 43 i 44 . By assuming that only one mode, e.g., the 
totally symmetric mode penetrates into the cavity for 
energy scales below A*, Eq. Q16p formally reproduces 
a one-channel Kondo model. Here, the goal of this Ap- 
pendix is to provide an alternative derivation, following 
the same terminology as the one used in Refs. 15814641 but 
keeping the channel-mixing terms. 

First, it is relevant to observe that the partition func- 
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tion in Eg . (|2"5)) can be expanded in powers of t as follows The notation (. . .}g implies that the the expectation 

value is taken with respect to the (free) actions of the 
electrons in the lead and dot. 



£ e ^nN J m [Vc] [VK] exp [_ 5o el _ gfi 

x ' / i Vm 



n— — oo 



m! 

m=0 



E E { -^r^ nNo J me- 3 " ((s int nt 



n——oo rn—0 



(Bl) 

Let us simplify the term ((S; n t) ) ■ 



{(S- mt ) m ) C a = t m I dn... dr m ^ ■■■ E ({cpiPidk 1 a 1 e~ l ' l>1 + h.c.) . . . (c Pml3m d kmam e 

fcipicti^i k m p m a m /3 m 



**■» +h.c.)). 

(B2) 



We alter the notation, such that the subscript 0fc = </>(rfc) 
is now taken to indicate the time label of the <j> operator. 
It is implicitly assumed that the <f> fields obey twisted 
boundary conditions. This expression has contributions 



only when m = 21 is even, and furthermore contains an 
equal number of c and c operators (and as a corollary d 
and J operators). The most general term looks like 



A =t 21 f ' dT il ...dT il dT jl ...dT jl J2 ••• E e- i ^+- + ^i-^--H) 

X ((cpixfti^i^! (dk^a^CpHPH ••■^ ! « 3! C P Ji ft ! ))' ( B3 ) 

I 



where and (ji, ■ • ■ ,j n ) are permutations 

of (1,...,2Z). Noting that the propagator of 
the electrons in the dot and electrons in the 
lead are independent of the channel index, i.e 
(dk iai (Ti)d k2a2 (T 2 )) = S aia2 (d kl (Ti)d k2 (T 2 )) and 
(cp^A-ri^p^)) = Sp 1 p a {c Pl (n)cp 2 ('ni)), and fur- 



thermore assuming a continuous spectrum (i.e. large 
lead and dot) such that k (d kl (Ti)d k2 (T 2 )') and 

J2 PlP2 (*-pi ( T i) c P2 ( T 2)) are proportional to each other, we 
obtain 



A = t 21 (MM) 1 / dr l ...<i7iae- <Wl+ - + * , -* , + 1 — -*« ) ^ ^ 
pi-pu fei...fc 2! 

)(d kl ...d kl d h+1 ...d k2l ). (B4) 

I 

In the last line we use the fact that for every complete each /3 we get a factor M. Since there are I such sums, 

contraction of the above correlators, the number of free we get an overall factor (MM) . 
indices for a and /? reduces from 21 to I. Summing over 

each a we get a factor M and similarly summing over Thus, the partition function reduces to 
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v-J — _~ Jo J<t> 



(11) 

1=0 v ' n=-oo 



x <^c Pl . . . c Pl c Pl+1 . . . c P2l ) (dki . . . dk l dk l+1 ■ ■ ■ dk 2l ) 

P1---P21 k 1 ...k 2 l 



(B5) 



The next step is now to use a "charge" representation of the partition function. 

We then transform the <p dependent part of the partition function (within square brackets) to the charge represen- 
tation as follows 



n— — oo 

oo 



2irnN 



4>(P)=4>a+2Trn 

[DQ] 



[Dtp] exp 



£ r / dr 



|Q 2 + iQ^} + i2irnQ(/3) - + ... + <&- - ... - 2/ ) 



0; - <p/ + l — — — <P2ZJ 

(B6) 

where we have introduced the auxiliary bosonic field Q(t), such that Q(J3) — Q(0). Defining the source term 



J(r;n,. . . ,t 2 j) = [<5(r - n) + . . . + 5(r - t ; ) - <5(r - r /+ i) - .. . - 5(r - t 2/ )] , 
we can rewrite Eq. (|B6j) as 



(B7) 



[X>Q] V] e «2jrn[Q(0)+iVo] e -B c / ^ drQ 2 
n— — oo 



0(^)=0o+2t™ 



[T>4>] exp 



^ |(j - J(t;ti, . . . ,r 2 ;)| 



P>Q] J2 5(Qa + N + n)e- Ec fo drQ-- 6 Q - J( r; n , . . . , r 2i ) 

71 — — OO 

n=-oo J-No-n 



n— — oo 



i 2Z 

Q(r) + 7V + n-5]e(r-r 2 )+ ^ 6 ( r " T >) 

i=l i=i+l 
( 2i 

Q( T )_ n + ^e( r _ r .)_ £ e(r-rO 

i=l t=J+l 



(B8) 



r 



where in the last step we used the transformation Q(t) — > 
N -Q(t). 

We now focus on the charge degeneracy point iVo = 1/2 
and assume that we are at very low temperatures (i.e. 
(3E C > 1 or more precisely /3A* > 1). This implies 
that we can focus only on the lowest two charge (energy) 
states. This mean trajectories of Q(t) are restricted to 
access these two charge states Q(t) = 0, 1. This pro- 
jection to the lowest energy states causes the time la- 
bels in Eq. (|B5I) to be nested, i.e. the trajectory of the 
charge (in imaginary time) consists of a succession of 



blips. Let us focus on one such trajectory is given by 
the nesting n > t; + i > t 2 > r/ +2 > ... > r 2 ;, and re- 
name T;_|_i — > t{,...,t 2 ; = t[. It is simple to see that 
permuting the time labels of the various nestings give 
equivalent contributions. There are l\ x l\ equivalent con- 
tributions, corresponding to permutations of the labels 
(ti, . . . , Ti) and (r{, . . . , r[ ). Furthermore, the contribu- 
tion from the trajectories 0^1— > ^ . . . ^ and 
1— > — > 1 — > 1 are identical. Thus in this regime 
we get 



2e ^ ^ 



;=o 



(I 



IV2 





11 x l\ l dn 



I) 



dr[ 



o 



dn 



ii 



dr[ 



x E ( c Pi(n)---Cp ! (r ; )cp ; (r 1 )...Cp ; (r;)^4 1 (T 1 )...4 i (r ; )4 i (r{)...4 ; (r 2; )^. (B9) 

pl...pjPj ...p[ k 1 ...kik' 1 ...k' l 
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We recall the identities 

Y] Gpi P 2( T i ~ Ta) = V ( c PiC P2 > = ^0 lim 
* * * — * — 

P1P2 P1P2 

and 



tt/0 



Vq- 



p->cc sin (7t//3(ti - T2)) n - 72 

7T//3 



£ G fclfc2 (n - r 2 ) = £ (d P J P2 ) = -1 lim - _ ^ 



1 



(iT I '8(Tl — To)) Ti — To ' 

where v>q and 1/1 denote the electron density of states in the lead and dot respectively. Thus we obtain 



be ; f Tl /" r '- 1 r i 

Z = 2e--^ V {v oVl t 2 UM) dn dr[... dn / dr{ 

;=0 Jo Jo Jo Jo 



1 



(BlOa) 



(BlOb) 



-1 2 



W-^)(i-Tft)...(7f-7^) 



(Bll) 



where (Pi, . . . , P^) is a permutation of (1, ... , /) and (— l) p is the corresponding sign (Note : interchanging a pair 
indices in (!,...,/) gives a factor of —1). Neglecting the constant factor of 2e ^ Eq. (|B11[) reduces to 



1=0 



Z = J2 {v QVl t 2 NM) 1 / dn dr[... 1 1 d Tl drl 



Jo 





rii<j in-Tj) 


\\\m<n ( T m T n)] 


nU (r[- Tll ) 







(B12) 



Renaming r{ — > T2, T2 — > T3, . . ., r ; ' — > T21 and introducing an ultraviolet cutoff l/r c = A*, we can rewrite this 
expression as 



Z = ^r(v QVl t 2 AfM) 1 



1=0 







■I 


lo T c 


JO T c 





T21 



■ exp 



2^(-l) < +'log 

i<j 



(B13) 



Recalling the partition function of the completely anisotropic single-channel Kondo model (in the language of the 
Coulomb gas expansion) 



_ x - / J\_V \ I an I - " «'2 



2n rfdn f Tl - T < dr 2 



'c Jo 



Tc dr 2n 



■ exp 



2^(-l)^log 



i<j 



'1 'j 



(B14) 



we show that it is identical to Eq. (|B13[) . where v is the density of states of the electrons of the Kondo model. The 
relationship between the one-channel Kondo theory and the parameters of the quantum RC circuit can be easily 
deduced 



= v Q vit 2 AfM. 



(B15) 



Appendix C: Perturbation theory near perfect 
transmission 

d + = d) + d 
(L = i{d^ - d) 
ip + (x) = il4(x) +ip s (x) 

4-{x)=i[tl>t(x)+Mx)], ( C1 ) 

which obey the commutation relations {ipc^c 1 } = 
and {d^jd^} — 2<5^<. Here, we have used the symbols 
It is convenient to introduce the Majorana fermions £ (£')=±. Using this representation and integrating out 

the charge held away from x = 0, the partition function of 
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the system is given as a coherent state functional integral We rewrite the action S = So + Sint using Eqs. (1551) . 

and (gBJ) as 

Z = J 'vfcVip+Thp-Vd+Vd-expl-S]. (C2) 



So — >Sq 



c=± 



fix 



and 



S int = / /> dTY;(|Ai| + C|Aa|)(5 c (r)Vc(0,r)(i_ c (T) := / drif int (r). 
Jo * ± Jo 



(C3) 



(C4) 



We define the constants = 2X^/vf, which describe the effective widths of the resonances of the (impurity) 
Majorana fermions. For the Majorana description near weak tansmission it is necessary to introduce the infrared 
cutoff r = max{r + ,r_}. This corresponds to the lifetime of the shorter lived Majorana fermion. The local (i.e at 
x = 0) Green's function of interest for the noninteracting Hamiltonian can be computed from the equations of motion 
on the lines of Ref. H3, and are given by 





rP 




/ dre 




Jo 




f P 


C-cK) = - 


/ dre 


Jo 




rP 


D ( (u> n ) = - 


/ dre 




Jo 



rJre^ r (r r [^(0,T)^(0,0)]) 



u> n sgn(w„) 



v F (oj n + rsgn(w„)) 



dTe^(r T [^-c(0,T)d c (0)]) 



i2X_ 



C 



sgn(w„) 



v F (w„ + Tsgn(w„)) 
2 



u) n + rsgn(w„) 



(C5) 
(C6) 
(C7) 



Note, w„ = (2n + 1)tt//3 denote (fermionic) Matsubara frequencies. Next, we compute the leading order correction to 
the charge response function 



K(t) = (T t [N(t)N(0)]) = ^(T t &(t)&(0) )=K (t) + K 1 (t) + K 2 (t) + O(X 3 ) 



(C8) 



where we have defined the functions 



K (t) = 



K x {r) = - 







dsi {T T 



^(r)^(0)ffi„t(ai) 



K 2 {t) = -^\ d Sl ds 2 (T T ^ c {T)^ c (0)H ia t(si)H int ( S2 ) 



(C9) 



Here, the subscript (. . .)o implies that the expectation value is taken with respect to the quadratic action 5*0. The 
contribution Kq{t) was computed in Eq. ([55)) from which we inferred R q = h/(2e 2 ). Let us explicitly denote the 
operators Cq as Cq = u^4> c + v c4> 2 c - The corrections to the Green's function due to backscattering from the function 
-Ki(t), is given by 



Ki{iu n ) 



2nE c 



J] n c {Nof (n +Cr 2 ) 2 log 
C=± V 



(T T [Vc(0,*i)d-c(C)(0)]) o 



us 

1 T 2£ c 



(CIO) 
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Here, we assume that we are in the Coulomb-blockaded regime and identify E c to be the ultraviolet cutoff . Wc 
define the functions k + (N q ) := cos(7rAT ) and k_(N ) := sin(7riVo)- The energy scale of the Fermi velocity, which is 
lost in the bosonization process where an infinite bandwidth is assumed, is reintroduced by fixing the cutoff 4 *^. The 
backscattering amplitudes are then determined to be 



A1.2 = ri 



2avpEcj 



(Cll) 



The contribution from the function K 2 (t) involves 3 types of diagrams 



1 fP 

K2(iu> n ) = — / dTdsids2e lU3nl 
n Jo 



E 

c=± 

^E 

c=± 



V27TO 

Ai + CA 2 



T T [0 c (r)^ c (O)(7 c (s 1 )C c ( S2 )]>(r T [V>c(O,s 1 )^c(O ) s 2 )])<T T [ < i_ c ( Sl )d_ f (s 2 )]) 



E 

Cl,2=± 



Ai + C1A2 \ / Ai + C2A2 



r T [^ c (T)^ c (o)c c (* 1 )Q(* a )])<r r [^c(o,«i)d_ c (« a )])<r r [^c(o,« 2 )d-c(*i)]) 



'2-ko, 



'2na 



T T [^(T)^ c (0)Q 1 (fli)C Ca ( S2 )]^<T r [^ 1 (0 ) «i)d- Cl (»!)]> (T T [V> C2 (0,s 2 )d_ C2 ( S2 )]> 

(C12) 



It is straightforward to show that the last term, which involves a product of cross-correlators of the impurity and lead 
Majorana fermions, is 0(A 4 ) and can be neglected. Furthermore, in this regime we should neglecting correlators in <j> c 
which involve more than a product of four <j) c fields, since the charge is pinned due to strong Coulomb interactions. 
Thus we obtain 



K 2 {iu n ) = - 



7 



2nE r 



E, 



21 



log ^ - 1 + - log 1 + ^ 



u;, 



(C13) 



In our calculations we make use of the Fourier transform of the <j> l c propagator 







F c {iv k ) = - I dTe WkT (T T ^(r)^(0) 



2\v k \+2E a /ir' 

where Vf. = 2kn//3 denotes (bosonic) Matsubara frequencies. Wick rotating icu n — > lu + iS we obtain 



(C14) 



K 2 (u) = 



7 



2nE c 

7 



x l 2E 



C=± 



^ Kc (iV ) 2 (r 1+ Cr 2 ) 2 log(^ 



2nE r 



2E C 



^ K _ c (7V ) 2 (r 1+ Cr 2 ) 5 



E, 



2r 



log l-r - 1 + log 



(C15) 
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